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- . Abstract 

o . 

7-H ■ A "microscopic" justification of the "symmetric damping" model of a quantum oscillator with time-dependent 

frequency and time-dependent damping is given. This model is used to predict results of experiments on 
simulating the dynamical Casimir effect in a cavity with a photo-excited semiconductor mirror. It is shown that 
the most general bilinear time-dependent coupling of a selected oscillator (field mode) to a bath of harmonic 
oscillators results in two equal friction coefficients for the both quadratures, provided all the coupling coefficients 
h— j , are proportional to a single arbitrary function of time whose duration is much shorter than the periods of all 

1 oscillators. The choice of coupling in the rotating wave approximation form leads to the "mimimum noise" 

' model of the quantum damped oscillator, introduced earlier in a pure phenomenological way. 

p_ _ n , Keywords: Dynamical Casimir effect, nonstationary quantum damped oscillator, Heiscnberg-Langevin equations, 
' noncommuting noise operators, ultrashort laser pulses, Wigner function, general bilinear bosonic coupling, rotating 
wave approximation. 
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1 Introduction 

Semiconductor mirrors have been used in the laser physics since the first years of the laser era [Uf4| . Recently 
they have found once more application. Namely, it is quite probable that they will help to realize a long-standing 
dream of an experimental verification of the so called Dynamical Casimir Effect (DCE). The story of this effect 
\ began 40 years ago, when Moore [5] showed that the motion of ideal mirrors could result in a creation of quanta 
of the electromagnetic field from the initial vacuum state. Since then, this exciting phenomenon was a subject of 
numerous theoretical studies (see reviews in [BHH]). 

The effect seemed to be extremely small, since the velocities of boundaries achievable in a laboratory are much 
less than the velocity of light. For this reason, it was thought for a long time that hardly the DCE could be 
observed. The situation began to change in 1990s, when a possibility of an enhancement of the effect inside cavities 
with oscillating boundaries under the conditions of parametric resonance was discovered [9TH3j. But exciting high 
frequency oscillations of rigid walls (about few GHz or higher) and maintaining them for a sufficiently long time 
(about 1 s) [10] is still a great experimental challenge. 

However, the task can be simplified, if one recognizes that all schemes that could be used to create quanta from 
the initial vacuum state of some circuit or cavity are based, after all, on a parametrical change of the resonance 
frequency of the system. But this can be achieved not only by changing the geometry (moving walls), but by 
changing other parameters, for example, the inductance and capacitance of some quantum circuit or properties of 
some Josephson junction |14||T5]. Nowadays this idea has been rediscovered and embodied in different practical 
schemes [TfJMT] . 

At the same time the idea of simulating "nonadiabatic Casimir effect" and other quantum phenomena using a 
medium with a rapidly decreasing in time refractive index ( "plasma window" ) was formulated by Yablonovitch [22) . 
Different theoretical schemes based on fast changes of the carrier concentration in semiconductors illuminated by 
laser pulses were discussed in [23H25] , and a possibility of creating an effective semiconductor microwave mirror in 
this way was confirmed experimentally |26j . The further development was a proposal (named "MIR experiment" ) 
[27j to simulate a motion of a boundary using an effective electron-hole "plasma mirror," created periodically on 
the surface of a semiconductor slab (attached to some part of the superconducting wall of a high-Q cavity) by 
illuminating it with a sequence of short laser pulses. The amplitude of an effective displacement of the boundary 
is determined in such a case by the thickness of the semiconductor slab, which can be made up to few millimeters, 
resulting in a relative change of the cavity eigenfrequency up to 10~ 2 , instead of the maximal possible value 10 -8 
for mechanically driven mirrors (this limit is due to tremendous internal stresses arising inside a material if the 
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frequency of the surface vibrations belongs to the GHz band [T0J). The current status of the MIR experiment 
was reported in [28] . A possibility to enhance and control the static Casimir force by illuminating semiconductor 
surfaces was shown in [29,50). 

Note that the thickness of the photo-excited conducting layer on the surface of the semiconductor slab is much 
smaller than the thickness of the slab itself. It is determined by the absorption coefficient of laser radiation, so it is 
about few micrometers or less, depending on the laser wavelength. Therefore laser pulses with the surface energy 
density about few /zJ/cm 2 can create a highly conducting layer with the carrier concentration exceeding 10 17 cm~ 3 , 
which gives rise to an almost maximal possible change of the cavity eigenfrequency for the given geometry. On 
the other hand, the conductivity of the layer is not extremely high due to a moderate value of the mobility in 
the available materials (such as highly doped GaAs), which is about lm 2 V s -1 [35]. For this reason, effects of 
dissipation during the excitation-recombination process inside the layer cannot be neglected, since they can change 
the picture drastically [7j. 

A model taking into account the dissipation was introduced in [71131] . It was based on the assumption that the 
effects of dissipation can be described within the frameworks of the linear Heisenberg-Langevin operator equations 
of the form (h = 1) 

dx/dt = p-j x (t)x + F x (t), dp/dt = - lp {t)p-uj 2 {t)x + F p {t), (1) 

where x and p are dimensionless quadrature operators of the selected mode of the EM field. These operators are 
normalized in such a way that the mean number of photons equals Af — \(p 2 + x 2 — 1). Two noise operators F x (t) 
and F p (t) with zero mean values (commuting with x and p) are necessary to preserve the canonical commutator 
[x(t),p(t)} — i. These operators were assumed to be delta-correlated, 

(F j (t)F k (t'))=S(t-t') Xjk (t), j,k = x,p. (2) 

The system of linear equations ([1]) can be solved explicitly for arbitrary time-dependent functions 7 XiP (i), ui(t) 
and F XiP (t). It appears that the commutation relation [x(t),p(t)\ = i can be preserved exactly under the condition 
Xxp — Xpx = 2i7(t) for any function j(t) = (j p + j x ) /2 [7J|3T]. The "symmetric damping" model with j x — ~f p 
was chosen in |JJ[3T] mainly due its simplicity. Later, the generic case of j x ^ j p was studied in [32 -34) . It was 
shown that there are some phenomenological reasons to connect the damping and noise coefficients as follows: 
Xxp = ~Xpx = ij(t)t Xx X = l' X G, Xpp = IpG, where the factor G > 1 depends on the reservoir temperature. 
However, the concrete value of the "asymmetry coefficient" y = (7^ — 7 X ) / (jp + 7^) remained as a free parameter. 
Although the choice y = permits one to simplify immensely many formulas |34j , it seems that there is no possibility 
to deduce the value of y from some general phenomenological principles [33 s , except for "aesthetic" ones. 

On the other hand, the physical results, such as the rate of photon generation, can depend significantly on the 
value of y 32] . This is because all formulas for the physical quantities contain a special solution of the classical 
nonstationary oscillator equation e-\-uj 2 j{t)e = with 0J 2 f{t) = w 2 (t) + S(t) — S 2 (t), where S(t) = [^/ x {t) — J P (t)} /2. 
If the damping coefficient "f(t) is much smaller than the amplitude of variations of the cavity eigenfrequency uj(t) 
and if this coefficient varies slowly with time (or it is constant), then the corrections to the effective frequency, 
5(t) and S 2 (t), can be safely neglected. This case was considered in [35], where it was shown explicitly that only 
the total damping coefficient 7 is important if "f x = const, j p = const and 7 <C uj. But in the case of the MIR 
experiment the total damping coefficient varies very fast, in the time scale of the order of the recombination time, 
which must be less than 30 ps [7J, whereas the period of the field mode is about 400 ps. Moreover, the maximum 
value of 7 in this case is only twice smaller than the amplitude of variations of u>(t). Therefore the term 5(t) cannot 
be neglected, and it gives rise to significant changes in comparison with the symmetric case of y — S(t) = [33)l34j . 

In principle, the value of y (as well as of all the damping and noise coefficients) could be deduced from some 
rigorous model, which would take into account explicitly (i) the coupling of the field mode with electron-hole 
pairs inside the semiconductor slab, (ii) the coupling of electrons and holes with phonons or other quasiparticles, 
responsible for the damping mechanisms, and (iii) a time dependence of the number of carriers, which disappear after 
a short recombination time. Unfortunately, it seems that no model of this kind is available now. Therefore I consider 
in this paper a surrogate "microscopic" model, where the real dissipativc system is replaced by a set of harmonic 
oscillators (bosonic reservoir) and the real interactions are replaced by an effective quadratic bilinear bosonic 
Hamiltonian. A justification is that this simple model was used in many text books when the authors considered 
field modes in non- ideal cavities [36H39] . Actually, it was considered in numerous publications for a long period since 
at least 1950s [40H5T] . The difference from the numerous existing studies is that here I consider the most general 
bilinear coupling with time- dependent coefficients. In the overwhelming majority of other papers only some specific 
couplings were considered, such as the so-called RWA (Rotating Wave Approximation) coupling 44 46 48,49 5 lj, the 
coordinate-coordinate coupling [40 42,45,47,51], the momentum-coordinate coupling 40, 5T], a sum of coordinate- 
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coordinate and momentum-coordinate couplings |50) . but all the coefficients in the interaction Hamiltonians were 
supposed to be time-independent. 

The scheme adopted in this paper was introduced in |52) and developed in |53) . It was used for studying the 
dccoherence and transfer of quantum states between modes in the Fabry-Perot cavity with resonantly oscillating 
boundaries in [54] . The next section is devoted to a generalization to the case of arbitrary time-dependent coupling 
coefficients. In Sec. [3] I show how the two equal damping coefficients j x = j p arise naturally in the special case 
of short-time interaction (which corresponds perfectly to the conditions of the MIR experiment), if all coupling 
coefficients are proportional to a single function of time (which can be connected with the total number of carriers 
inside the photo-excited semiconductor film). The "minimum noise" model arises naturally under the additional 
assumption of the RWA coupling. 



2 Oscillator coupled to nonstationary bosonic bath 



Let us assume that the "central" oscillator with the time-dependent frequency to(t), described in terms of the 
"coordinate" xq and canonically conjugated "momentum" po (representing the selected field mode in the cavity) is 
coupled to a large number of the "bath" oscillators with constant frequencies u>i, and coordinates Xi and momenta 
Pi. The Hamiltonian of the whole (closed) system has the form 



1 1 N N 

H = -[pl+ uj 2 (t)xl] +r^(p, 2 + <Axf) + 2J (zipipo + VipiXQ + UiXipo + giXtxo) , 



(3) 



where coupling coefficients Zi, Vt, u i: gi can be arbitrary functions of time. The effective masses of every oscil- 
lator are assumed to be equal to unity (this can be achieved by rescaling the coordinates). It is convenient to 
introduce the (2N + 2)-dimensional vector q = (Q,£), whose components are defined as follows: Q = (po,Xq), 

£ = (pi,P2, ■ ■ ■ iPN, Xl,X2, ■ ■ ■ ,X N ). 

The simplest description of the evolution of quantum systems with multidimensional quadratic Hamiltonians, 
such as ([3]), can be achieved in the Wigner representation |55 [ I56 ) . The time-dependent Wigner function of the 
whole system can be written as 

(4) 



W(q,t) = J G(q,q / ,t)W(q',0)dq / 



where the propagator G(q, q',£) satisfies the initial condition G(q, q',0) = <5(q — q'). One of many remarkable 
properties of the Wigner function is that its evolution in the case of quadratic Hamiltonians is governed by the first- 
order partial differential equation [56 ( 157 ] (whereas the Schrodinger equation or the equation for the density matrix 
contain the second order derivatives). Therefore, the whole propagator is reduced to a delta- function [57H61] . 
G(q, q',i) = 5(q — q*(i;q')), where vector q*(t;q') = R(i)q' is the solution to the classical equations of motion 
satisfying the initial condition q*(0; q') = q'. Recall that for systems with quadratic Hamiltonians, the Heisenberg 
equations of motion for operator q have exactly the same form as the equations for classical generalized coordinates 
(or the first-order mean values), so R(i) is a symplectic (2N + 2) x (2N + 2) matrix satisfying the linear equation 



R = ^R, 



R(0) 



12JV+2) 
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where 1m is the M x M unity matrix. It is useful to split the matrices R and A into rectangular blocks in 
accordance with the structure of vector q = (Q, £): 



R 



Rn 
R21 



R12 
R22 



A = 



An Aa 
A21 A22 



Then Eq. ([5]) is equivalent to the following set of equations: 



Rn 


— .A11R11 - 


I- .4i2R21, 


Rn(0) 


= 1, 


R21 


= .A21RU " 


- -422R21, 


R 2 i(0) 


-0, 


R12 


= >4iiRi2 - 


- -4l2R22, 


Ria(O) 


= 0, 


R22 


= ^2lRl2 " 


- ^22R22, 


R 22 (0) 


= I. 



For Hamiltonian (j3)), the matrices Ajk have the following explicit forms: 



An 







A 



22 





diag(l,...,l,...) 



-diag (ujf, . 







(6) 
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(8) 
(9) 
(10) 

(11) 
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. .) means the diagonal matrix with the elements ai, 


, . . , ttjj . 





(12) 



The notation diag (ai 

I assume that the initial total Wigner function is factorized with respect to the Q and £ variables and that the 
initial Wigner function of the reservoir is Gaussian (hereafter h = 1) 



W(q,0) = W (Q)Wi(0, Wx(0 = (detF)- 1 / 2 exp 



1 



-.N 



1. 



(13) 



Here F is a symmetric 2N x 2N positively definite matrix. The first-order mean values of the reservoir quadratures 
are supposed to be equal to zero, for the sake of simplicity. In particular, the state (fT3")) may correspond to some 
mixed thermal state or some pure squeezed vacuum state. The dynamics of the subsystem Q is described by the 
reduced Wigner function W(Q,t) = J W(Q, f, t)d 2N £/(2tt) n , which can be written as follows, 



w(Q,t) = (2tt)- n J dQ'de'deWo(Q')WiK')*(Q - R nQ' - RiaO*(£ - R 2iQ' - R 22^')- 



(14) 



The integration over c?£ simply removes the second delta-function. To perform the integration over d£' one should 
replace the first delta-function with its integral representation <5(x) = J e lkx d 2Ar k/(27r) 2Ar . Thus one arrives at two 
Gaussian integrals (the first one over d£' and the second over dk), which can be easily calculated exactly. The 
final result is the formula, which has the same form as (2]), but with the variables Q and Q' instead of q and q', 
and with the averaged propagator 



G(Q,Q',t) = (27r)- n [dctM4<)]- 1 / 2 exp 
where symmetric 2x2 matrix Al*(t) equals 

M*{t) =Ria(t)FRi 2 (t) 



--(Q-RnQOM^^Q-RiiQ') 



(15) 



(16) 



A direct inspection shows that the propagator (|15p (consequently, its convolution with any initial function) 
satisfies the Fokker-Plank equation 



dW d d 2 W 

at oQ a H oQadQp 

with the following time-dependent drift matrix A and time-dependent matrix of diffusion coefficients D: 

A = RiiR^ = An + A12R21R11 , 
2D = M* - AM, - M,A = A 12 (R22 - RaiR^Ria) FR12 + R12F (r 22 - RaiRn Ria) M2 



(17) 

(18) 
(19) 

On the other hand, it is well known 56 ,62 ,63] that the description of the evolution of an open quatum system 
by means of the equation (JTTJ) is equivalent to the Heisenberg-Langevin equation 



dQ/dt = AQ + x(t) 
with delta-correlated in time noise operators Xa(t), 

(x a (t)) = 0, (x a (t)Mt')) =s(t- x aP . 

The 2x2 noise matrix X and its transposition X are related to the diffusion matrix D as follows, 



4D = X + X, X = 



Xn X12 
X21 X22 



(20) 
(21) 

(22) 
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In many applications (including the case considered here) the elements of the interaction matrices A12 and A21 
are small. Then one can use the method of successive perturbations. In the zeroth approximation the solution of 
Eq. (fTOf reads (as soon as matrix A22 is supposed to be time independent) 



B^(t) = exp{A 22 t) 



diag(cos u>it) diag(—Wi sin Wji) 
diag(wr 1 sincjjt) diag(coswji) 



(23) 



Putting H 2 °2 (t) into the right-hand side of ([5J one obtains the first-order solution for the matrix R2 



R$ (t) = exp (A 2 2t) I exp (-A 22 t) A 21 (t)K^ (r) dr. (24) 
Jo 

Within the same accuracy, formula (fT9]l can be simplified as 

2D = .A12R22 FR^ + R^FR^Uu. (25) 

3 Special case of very short interaction time 

At this point it is necessary to take into account one of the most important results of the previous theoretical 
studies on the DCE in cavities with semiconductor mirrors, namely, that the recombination time of the carriers 
must be much shorter than the period of the field oscillations of the fundamental cavity mode and the duration 
of each laser pulse must be of the order of the recombination time or shorter. For the frequency luq « 2.5 GHz 
the period of oscillations is about 400 ps, while the recombination time should not exceed 35 ps [7j- This means 
that during the whole interval of time, when u>(t) is (slightly) different from the initial (and final) value coq and 
the coupling coefficients arc different from zero, the matrix (£) = exp (.A22O m Eq. (f2~4"|) . as well as the matrix 
Rix (t) (although it cannot be written as a simple matrix exponential function if the frequency w depends on time) , 
can be replaced by the unity matrices. Of course, one has to make an additional (but quite realistic) assumption 
that the coupling coefficients with the phonon modes with co* 3> are negligibly small. 

Then Eq. Q18|) leads to the following first-order approximation for the drift matrix A governing the evolution 
of the average values of the subsystem variables: 



A(D (t) = A n (t) + A 12 (t)R 2 1 1 ) (t) = An (t) + A 12 (t) /* A 21 (r) dr. 

Jo 



(26) 



Similarly, one obtains the approximate formula R^ (t) ~ J Q * Ai 2 (t) dr, so that the diffusion matrix (|25[) can be 
written as 

2Ty = A 1 2(t)FR[ 1 l(t)+R[ 1 2 \t)FA 1 2(t)=A 1 2(t)F i 12 (r)dr+/ A 12 (r) drF A 12 (t). (27) 

Jo Jo 

Taking into account the explicit form (fT2|) of matrices Ai 2 (t) and A 2 i(t), the elements of the dissipative part 
of the drift matrix ^ 

fx(t) = A (1) (t) - An(t) = Au{t) f A 2 i{t) dr (28) 

Jo 

can be written as follows, 



N .t N t 

Mil(*)=y^/ dr [v k (t)u k {T) - g k (t)z k (T)] , fj, 12 (t) = V / dr [v k (t)g k {T) - g k (t)v k (T)} , 



^2i{t) I dr [u k {t)z k (T) - z k (t)u k (T)] , 



(*) = V" / dr [u k (t)v k (T) - z k (t)g k (r)} 

L.-1 JO 



fc=l" U fe=l 

In the case of the DCE with a semiconductor mirror, the time dependence of the coupling coefficients between 
the field and the "reservoir" is introduced in order to take into account the time dependence of the number of 
carriers inside the semiconductor (since there is no coupling in the absense of carriers, if the quality factor of the 
empty cavity is big enough). Therefore I assume that all coupling coefficients are proportional to the single time 
dependent factor v(t) (related somehow to the number of carriers): 

u k {t) = v[t)U k , v k (t) = v(t)V k , g k (t) = v(t)G k , z k (t) = v(t)Z k , (29) 
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where the coefficients Uk,Vk,Gk, Zk do not depend on time. Then straightforward simple calculations lead to the 
following remarkable result: 



N „t 
M12 = M21 = 0, nn = [i 2 -2 = A(t) y2 {UkVk - G k Z k ) , A(t) = v(t) / v{r)dT. 

k=i J ° 



(30) 



Consequently, two equal damping coefficients in Eq. ([1} arise automatically: "f x = /122 = Mil = 7j>- 

diag^f/i) 

diag (/,; 

an equilibrium variance matrix for the reservoir variables), then Eqs. H2), (ITT)) and ([25]) result in the following 
expressions for the elements of the 2x2 symmetric diffusion matrix D: 



If the reservoir variance matrix has the diagonal form F 



(in particular, F may be 



N 



D u (t) = \(t) fk HV k 2 + Gl) , D 22 (t) = X(t) fk HZ 2 k + U 2 k ) , 



k=l 



k=l 



N 



D 12 (t) = -X(t) fk HV k Z k + G k U k ) 



k=l 



These coefficients become especially simple if the interaction Hamiltonian in §5§ has the RWA form: 



N 

Hint = v(t) ^2 (pkaoal 



p* k a k al 



ujqXq + ip A Lo k Xk + ipk 
ao = 7== — , a k = 



k=l 



'2uj 



'2uk 



where u)$ is the initial frequency of the central oscillator. Then 

Re(p fc ; 



G k = ^/uj a uj k Re(pk), Z k 



y/U UJk 



Vfc = ,— lm(pk), U k = -\— Im(p fc ), 
V w fc V UJQ 



(31) 
(32) 

(33) 
(34) 



so that D12 = and flu = = X(t)wo Ylk=i fk^k\pk\ 2 ■ Using (f22|) one arrives at the set of the noise 

coefficients which has exactly the same structure as the minimum noise set of coefficients introduced in [7|I31| in 
the framework of a phenomenological approach: 



Xxp 



■Xpx = ij(t), Xpp = ulxxx = l{t)u G. 



(35) 



The same set (|35() with 7 = const arises in the case of the RWA coupling with time-independent coefficients for 
w t 1 [S3]. This observation makes the "bridge" between the short-time and long-time limits. It supports the 
assumption made in [71131] that the coefficient G can be identified with the mean photon number in the equilibrium 
state with the given temperature T and frequency uiq, G — coth [hu /(2kBT)}. 



4 Conclusion 

It was shown that the "symmetric minimum noise" set of coefficients in the Heisenberg-Langevin equations for a 
quantum damped nonstationary oscillator, used in the studies of the dynamical Casimir effect in a cavity with a 
photo-excited semiconductor mirror, arises quite naturally from any bilinear interaction between the selected field 
mode and the bosonic bath. The main conditions are: (i) the duration of interaction is much shorter than the 
period of the oscillator and (ii) all the coupling coefficients are proportional to a single function of time. Of course, 
real interactions are much more complicated, so the construction of a more realistic model remains to be a very 
interesting and important problem. Nonetheless, the results obtained give a strong support to the conclusions of 
recent studies [Tl fSlIM] concerning the feasibility of the MIR experiment and the necessary conditions for its success. 



Acknowledgements 

A partial support from CNPq (DF, Brazil) is acknowledged. 



G 



References 

[1] P. G. Eliseev and M. A. Man'ko, "The use of a semiconductor mirror for Q-switching laser," Zh. Tekh. Phys., 
36, 2215-2216 (1966) [Sov. Phys. Tech. Phys., 11, 1657-1658 (1967)]. 

[2] C. H. Carmichael and G. N. Simpson, "Generation of giant optical maser pulses using semiconductor mirror," 
Nature, 202, 787-788 (1964). 

[3] W. R. Sooy, M. Geller, and D. P. Bortfeld, "Switching of semiconductor reflectivity by a giant pulse laser," 
Appl. Phys. Lett., 5, 54-56 (1965). 

[4] M. Birnbaum, "Modulation of the reflectivity of semiconductors," J. Appl. Phys., 36, 657-658 (1965). 

[5] G. T. Moore, "Quantum theory of the electromagnetic held in a variable-length one-dimensional cavity," J. 
Math. Phys., 11, 2679-2691 (1970). 

[6] V. V. Dodonov, "Nonstationary Casimir effect and analytical solutions for quantum fields in cavities with 
moving boundaries," in: M. W. Evans (Ed.), Modern Nonlinear Optics, Advances in Chemical Physics Series, 
Wiley, New York (2001), Vol. 119, Pt. 1, pp. 309-394 [quant-ph/010608l] . 

[7] V. V. Dodonov and A. V. Dodonov, "Quantum harmonic oscillator and nonstationary Casimir effect," J. Russ. 
Laser Res., 26, 445-483 (2005). 

[8] V. V. Dodonov, "Dynamical Casimir effect: some theoretical aspects," J. Phys.: Conf. Ser., 161, 012027 
(2009). 

[9] V. V. Dodonov, A. B. Klimov, and V. I. Man'ko, "Generation of squeezed states in a resonator with a moving 
wall," Phys. Lett. A, 149, 225-228 (1990). 

[10] V. V. Dodonov and A. B. Klimov, "Generation and detection of photons in a cavity with a resonantly oscillating 
boundary," Phys. Rev. A, 53, 2664-2682 (1996). 

[11] A. Lambrecht, M.-T. Jaekel, and S. Reynaud, "Motion induced radiation from a vibrating cavity," Phys. Rev. 
Lett, 77, 615-618 (1996). 

[12] D. A. R. Dalvit and F. D. Mazzitelli, "Creation of photons in an oscillating cavity with two moving mirrors," 
Phys. Rev. A, 59, 3049-3059 (1999). 

[13] G. Plunien, R. Schiitzhold, and G. Soff, "Dynamical Casimir effect at hnite temperature," Phys. Rev. Lett., 
84, 1882-1885 (2000). 

[14] V. V. Dodonov, O. V Man'ko, and V. I. Man'ko, "Correlated states in quantum electronics (resonant circuit)," 
J. Sov. Laser Res., 10, 413-420 (1989). 

[15] V. I. Man'ko, "The Casimir effect and quantum vacuum generator," J. Sov. Laser Res., 12, 383-385 (1991). 

[16] S. De Liberato, C. Ciuti, and I. Carusotto, "Quantum vacuum radiation spectra from a semiconductor micro- 
cavity with a time-modulated vacuum Rabi frequency," Phys. Rev. Lett., 98, 103602 (2007). 

[17] K. Takashima, N. Hatakenaka, S. Kurihara, and A. Zeilinger, "Nonstationary boundary effect for a quantum 
flux in superconducting nanocircuits," J. Phys. A, 41, 164036 (2008). 

[18] K. Takashima, S. Matsuo, T. Fujii, N. Hatakenaka, S. Kurihara, and A. Zeilinger, "Squeezing of a quantum 
flux in a double rf-SQUID system," J. Phys.: Conf. Ser., 150, 052260 (2009). 

[19] A. V. Dodonov, "Photon creation from vacuum and interactions engineering in nonstationary circuit QED," 
J. Phys.: Conf. Ser., 161, 012029 (2009). 

[20] J. R. Johansson, G. Johansson, C. M. Wilson, and F. Nori, "Dynamical Casimir effect in a superconducting 
coplanar waveguide," Phys. Rev. Lett, 103, 147003 (2009). 

[21] S. De Liberato, D. Gerace, I. Carusotto, and C. Ciuti, "Extracavity quantum vacuum radiation from a single 
qubit," Phys. Rev. A, 80, 053810 (2009). 



7 



[22] E. Yablonovitch, "Accelerating reference frame for electromagnetic waves in a rapidly growing plasma: Unruh- 
Davies-Fulling-De Witt radiation and the nonadiabatic Casimir effect," Phys. Rev. Lett., 62, 1742-1745 
(1989). 

[23] E. Yablonovitch, J. P. Heritage, D. E. Aspnes, and Y. Yafet, "Virtual photoconductivity," Phys. Rev. Lett., 
63, 976-979 (1989). 

[24] T. Okushima and A. Shimizu, "Photon emission from a false vacuum of semiconductors," Jpn. J. Appl. Phys., 
34, 4508-4510 (1995). 

[25] Y. E. Lozovik, V. G. Tsvetus, and E. A. Vinogradov, "Parametric excitation of vacuum by use of femtosecond 
laser pulses," Phys. Scr., 52, 184-190 (1995). 

[26] C. Braggio, G. Bressi, G. Carugno, A. Lombardi, A. Palmieri, G. Ruoso, and D. Zanello, "Semiconductor 
microwave mirror for a measurement of the dynamical Casimir effect," Rev. Sci. Lustrum., 75, 4967-4970 
(2004). 

[27] C. Braggio, G. Bressi, G. Carugno, C. Del Noce, G. Galeazzi, A. Lombardi, A. Palmieri, G. Ruoso, and D. 
Zanello, "A novel experimental approach for the detection of the dynamic Casimir effect," Europhys. Lett., 
70, 754-760 (2005). 

[28] A. Agnesi, C. Braggio, G. Bressi, G. Carugno, F. Delia Valle, G. Galeazzi, G. Messineo, F. Pirzio, G. Reali, 
G. Ruoso, D. Scarpa, and D. Zanello, "MIR: an experiment for the measurement of the dynamical Casimir 
effect," J. Phys.: Conf. Ser., 161, 012028 (2009). 

[29] N. Inui, "Numerical study of enhancement of the Casimir force between silicon membranes by irradiation with 
UV laser," J. Phys. Soc. Jpn., 73, 332-339 (2004). 

[30] F. Chen, G. L. Klimchitskaya, V. M. Mostepanenko, and U. Mohideen, "Control of the Casimir force by the 
modification of dielectric properties with light," Phys. Rev. B, 76, 035338 (2007). 

[31] V. V. Dodonov, "Time-dependent quantum damped oscillator with 'minimal noise': application to the non- 
stationary Casimir effect in nonideal cavities," J. Opt. B, 7, S445-S451 (2005). 

[32] V. V. Dodonov and A. V. Dodonov, "The Heisenberg-Langevin model of a quantum damped harmonic oscil- 
lator with time-dependent frequency and damping coefficients," J. Russ. Laser Res., 27, 379-388 (2006). 

[33] V. V. Dodonov and A. V. Dodonov, "Theory of the dynamical Casimir effect in nonideal cavities with time- 
dependent parameters," J. Phys.: Conf. Ser., 99, 012006 (2008). 

[34] V. V. Dodonov, "Photon distribution in the dynamical Casimir effect with an account of dissipation," Phys. 
Rev. A, 80, 023814 (2009). 

[35] V. V. Dodonov, "Dynamical Casimir effect in a nondegenerate cavity with losses and detuning," Phys. Rev. 
A, 58, 4147-4152 (1998). 

[36] W. H. Louisell, Quantum Statistical Properties of Radiation, Wiley, New York (1973). 

[37] M. O. Scully and M. S. Zubairy, Quantum Optics, Cambridge University Press, Cambridge (1997). 

[38] U. Weiss, Quantum Dissipative Systems, World Scientific, Singapore (1999). 

[39] C. W. Gardiner and P. Zoller, "Quantum Noise," Springer, Berlin (2000). 

[40] M. Toda, "On the theory of the Brownian motion," J. Phys. Soc. Jpn., 13, 1266-1280 (1958). 

[41] T. Kotera and M. Toda, "On the theory of Brownian motion and spin relaxation," J. Phys. Soc. Jpn., 14, 
1475-1490 (1959). 

[42] V. B. Magalinskii, "Dynamical model in the theory of the Brownian motion," Zh. Eksp. Teor. Fiz., 36, 
1942-1943 (1959) [Sov. Phys. - JETP, 9, 1381-1382 (1959)]. 

[43] C. George, "Mouvement Brownien d'un oscillateur quantique," Physica, 26, 453-477 (1960). 

[44] J. P. Gordon, L. R. Walker, and W. H. Louisell, "Quantum statistics of masers and attenuators," Phys. Rev., 
130, 806-812 (1963). 



8 



[45] W. G. Ford, M. Kac, and P. Mazur, "Statistical mechanics of cnscmblics of coupled oscillators," J. Math. 
Phys., 6, 504-515 (1965). 

[46] A. E. Glassgold and D. Holliday, "Quantum statistical dynamics of laser amplifiers," Phys. Rev., 139, A1717- 
A1734 (1965). 

[47] P. Ullersma, "An exactly solvable model for Brownian motion. I. Derivation of the Langevin equation," Physica, 
32, 27-55 (1966). 

[48] B. R. Mollow, "Quantum theory of field attenuation," Phys. Rev., 168, 1896-1919 (1968). 

[49] R. J. Glauber, "Coherence and quantum detection," in: R. J. Glauber (Ed.), Quantum Optics, Proceedings 
of the International School of Physics "Enrico Fermi," Course XLII, Varenna, 1967, Academic, New York 
(1969), pp. 15-56. 

[50] G. S. Agarwal, "Brownian motion of a quantum oscillator," Phys. Rev. A, 4, 739-747 (1971). 

[51] G. W. Ford, J. T. Lewis, and R. F. O'Connell, "Quantum Langevin equation," Phys. Rev. A, 37, 4419-4428 
(1988). 

[52] V. V. Dodonov, "Different approaches to the problem of dissipation in quantum mechanics," J. Kor. Phys. 
Soc, 24, S70-S74 (1991). 

[53] V. V. Dodonov, O. V. Man'ko, and V.I. Man'ko, "Quantum nonstationary oscillator: models and applications," 
J. Russ. Laser Res., 16, 1-56 (1995). 

[54] V. V. Dodonov, M. A. Andreata, and S. S. Mizrahi, "Decoherence and transfer of quantum states of field 
modes in a one-dimensional cavity with an oscillating boundary," J. Opt. B, 7, S468-S479 (2005). 

[55] M. Hillery, R. F. OConnell, M. O. Scully, and E. P. Wigner, "Distribution functions in physics: fundamentals," 
Phys. Rep., 106, 121-167 (1984). 

[56] V. V. Dodonov, "Parametric excitation and generation of nonclassical states in linear media," in: V. V. 
Dodonov and V. I. Man'ko (Eds.), Theory of Nonclassical States of Light, Taylor & Francis, London (2003), 
pp. 153-218. 

[57] J. E. Moyal, "Quantum mechanics as a statistical theory," Proc. Cambridge Phil. Soc, 45, 99-124 (1949). 

[58] T. Takabayasi, "The formulation of quantum mechanics in terms of ensemble in phase space," Prog. Theor. 
Phys., 11, 341-373 (1954). 

[59] J. G. Kriiger and A. Poffyn, "Quantum mechanics in phase space III: linear transformations," Physica A, 91, 
99-112 (1978). 

[60] V. V. Dodonov, V. I. Man'ko, and V. N. Rudenko, "Quantum properties of high-Q macroscopic resonators," 
Sov. J. Quantum Electron. 10 1232-8 (1980). 

[61] M. Gadclla, J. M. Gracia-Bondia, L. M. Nicto, and J. C. Varilly, "Quadratic Hamiltonians in phase-space 
quantum mechanics," J. Phys. A, 22, 2709-38 (1989). 

[62] M. Lax, "Quantum noise. IV. Quantum theory of noise sources," Phys. Rev., 145, 110-129 (1966). 

[63] H. Haken, "Cooperative phenomena in systems far from thermal equilibrium," Rev. Mod. Phys., 47, 67-121 
(1975). 



9 



